Levy— Brownian motion on finite intervals: Mean first passage time analysis. 
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We present the analysis of the first passage time problem on a finite interval for the generalized 
Wiener process that is driven by Levy stable noises. The complexity of the first passage time 
statistics (mean first passage time, cumulative first passage time distribution) is elucidated together 
with a discussion of the proper setup of corresponding boundary conditions that correctly yield 
the statistics of first passages for these non-Gaussian noises. The validity of the method is tested 
numerically and compared against analytical formulae when the stability index a approaches 2, 
recovering in this limit the standard results for the Fokker-Planck dynamics driven by Gaussian 
white noise. 
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I. INTRODUCTION 

Stochastic Levy processes serve as paradigms for the 
description of many unusual transport processes leading 
to anomalous diffusion as characterized by an anomalous 
mean squared displacement, i.e., 



((x(t)-(x(t))) 2 ) = ((Ax) 2 )^t^ 



(1) 



which deviates with v ^ 1 from the linear dependence 
((Ax) 2 ) oc t that characterizes normal diffusion. In the 
above formula v stands for the anomalous diffusion ex- 
ponent that specifies the process at hand as either be- 
having subdiffusive (with < v < 1) .suverdiffusive 
(1 < v) or ballistic (for v = 2) 0, H 0, ||. Among 
the class of Levy processes, the free Levy flights (LF) 
represent a special class of discontinuous Markovian pro- 
cesses, for which the mean squared displacement, as de- 
fined in equation (JIJ diverges due to the heavy-tail dis- 
tribution of the independent increments Ax. In this case, 
the mean squared displacement is always superdiffusive; 
i.e. ((x(t) - (x(t))) 2 ) oc t 2 / a , where < a < 2 denotes 
the stability index of the Levy-Brownian motion process, 
see below in Sect. II. In contrast to the spatio-temporal 
coupling characterizing general forms of non-Mar kovian, 
or more precisely, semi-Markovian Levy walks (LW) 
[B S H Levy flights correspond to the class of Markov 
processes that emerge from a Langevin equation with 5- 
correlated, white Levy noise. Because LFs typically pos- 
sess a broad jump length distribution with an asymptotic 
power law behavior their trajectories display at all scales 
self-similar clustering of local sojourns that become inter- 
rupted by long jumps into the other location in the phase 
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space where a new clustering forms. Early discoveries of 
LF-like phenomena were related to intermittent chaotic 
systems and description of the motion of the fluid par- 
ticles in fully developed turbulence 0,00 Nowadays, 
their applications range from description of the dynam- 
ics in plasmas, diffusion in the energy space, self-diffusion 
in micelle systems and transport in polymer systems un- 
der conformational motion [J, [5j to the spectral analysis 
of paleoclimatic [f| or economic data @. Despite the 
ubiquitous use of LFs as phenomenological models for 
noise sources, their influence on the kinetics subjected to 
boundary conditions has been addressed scarcely only. 
For free normal diffusion, the knowledge of the Green 
function, together with the local boundary conditions, is 
sufficient to determine the first passage time statistics. 
The same information can be also obtained by use of the 
method of images or by solving the corresponding, local 
boundary value problem of the diffusive Fokker-Planck 
equation H,0|. 



In particular, for processes driven by white Gaussian 
noise these boundary conditions for reflection or absorp- 
tion are locally defined and well known |9lll0l|. However, 
for the case of LFs, the method of images fails, yield- 
ing results that contradict the Sparre- Andersen theorem 
[Hill. By virtue of the latter [l^, for any discrete 
time random walk starting out at Xq ^ with the step 
length sampled from a continuous, symmetric distribu- 
tion, the first passage time density decays asymptotically 
as i~ 3 / 2 . In order to further explore the intricate prob- 
lem of proper boundary conditions for absorption and 
reflection for non-Gaussian white noise composed of in- 
dependent Levy flight increments, we elucidate with this 
work the situation of a free, overdamped Levy-Brownian 
motion, being restricted to a bounded domain of attrac- 
tion. 
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II. RESTRICTED LEVY— BROWNIAN MOTION 

We consider the free Brownian motion on a restricted, 
finite interval that is driven by Levy stable noise. More 
specifically, the dynamical evolution of a stochastic state 
variable x(t) is described in terms of the Langevin equa- 
tion 



dx(t) 
dt 



at) 



(2) 



where ((t) denotes a Levy stable white noise process 
which is composed of independent differential increments 
that are distributed according to the stable density with 
the index a; i.e. L a ^(£; a, \x = 0). Put differently, 
stands for the generalized white noise process which is 
obtained from the time derivative of the corresponding 
Levy-Brownian (Markovian) stable process. The param- 
eter choice a = 2 yields the usual, ^-correlated Gaus- 
sian white noise. In contrast, with < a < 2 the 
corresponding Levy white noise is generated from Levy- 
Brownian motion possessing discontinuous sample paths 
with infinite variance and the higher cumulants. Its stati- 
tiscal properties can be however characterized by frac- 
tional moments of order v which exist and are finte for 
v < a < 2 [l4|. Because it is composed of indepen- 
dent increments, this Levy stable white noise also con- 
stitutes a singular white noise process whose autocor- 
relation again is formally <5-correlated |l4|. Here, the 
parameter a denotes the stability index, yielding the 
asymptotic power law for the jump length distribution 
being proportional to The parameter a char- 

acterizes a scale, (3 defines an asymmetry (skewness) of 
the distribution, whereas \i denotes the location param- 
eter. We deal only with strictly stable distributions not 
exhibiting a drift regime; this implies a vanishing loca- 
tion parameter [i = throughout the remaining part 
of this work. For a =/= 1, the characteristic function 
4>{k) — e~ tk ^ L aj p(£; a,[i = 0)d( of an a-stable ran- 
dom variable £ can be represented by 



<j)(h) = exp -a a \k\ a (l - i/3sign(*) 



while for a = 1 this expression reads 



tan ■ 



4>(k) = exp 



-o-lfcl 1 



i(3— sign(fc) In \k 

7T 



(3) 



(4) 



The three remaining parameters vary within the regimes 
a e (0,2], (3 e [-1,1], cr G (0,oo). 

The stochastic differential equation in © yields nor- 
mal, free Brownian motion when a = 2, and free super- 
diffusion when a £ (0,2). The numerical integration of 
equation (J2J has been performed by use of standard tech- 
niques of integration of stochastic differential equation 
with respect to the Levy stable measures and studied 
by use of the Monte Carlo methods 0]. In particu- 
lar, the error bars visible in the figures were calculated 



using the bootstrap method. The position of a Levy- 
Brownian particle has been obtained by a direct integra- 
tion of Eq ( Et leading to the following approximation 



it) = / C(s)ds 



N-1 

£(a*) 1/b c< 



(•5) 



where Q are independent random variables dis- 
tributed with the probability density function (PDF) 
L a ,/3(C', c, [i = 0) and NAs = t. All our illustrations in 
this work are based on trajectory calculations sampled 
from the Langevin equation J3J). Absorbing boundary 
conditions have been realized by stopping the trajectory 
whenever it reached the boundary, or, more typically, it 
has jumped beyond that boundary location. The condi- 
tion of reflection has been assured by wrapping the hit- 
ting (or crossing) trajectory around the boundary loca- 
tion, while preserving its assigned length. The Appendix 
A provides some further details on the numerical scheme 
for stochastic differential equations driven by Levy white 
noise. We also like to emphasize that our results omit 
cases when a = 1 with (3 ^ 0. In fact, this parameter set 
is known to induce instabilities in the numerical evalua- 
tion of corresponding trajectories [l6l IrH ITsl IT^. I20I l2l| . 
Thus, for our numerics with ^ 0, the parameter value 
a = 1 has been excluded from the consideration. 



III. MEAN FIRST PASSAGE TIMES FOR 
a-STABLE NOISES 

Our main objective is the investigation of the mean 
first passage time (MFPT) of Levy-Brownian motion on 
a finite interval as illustrated with Fig.^ The boundaries 
at both ends will be assumed to be of either absorbing or 
also of the reflecting type. With generally non-Gaussian 
white noise the knowledge of the boundary location alone 
cannot specify in full the corresponding boundary condi- 
tions for absorption or reflection, respectively. In par- 
ticular, the trajectories driven by non-Gaussian white 
noise depict discontinuous jumps, cf. Figs. [3 and 0] be- 
low. As a consequence, the location of the boundary itself 
is not hit by the majority of discontinuous sample tra- 
jectories. This implies that regimes beyond the location 
of the boundaries must be properly accounted for when 
setting up the boundary conditions. Most importantly, 
returns (i.e. so termed re-crossings of the boundary lo- 
cation) from excursions beyond the specified state space 
back into this very finite interval where the process pro- 
ceeds must be excluded. Thus, the problem of proper 
formulation of boundary conditions in such cases poses 
an open challenge that has not been addressed with suf- 
ficient care in the prior literature |22L |24| . In contrast to 
the case with normal diffusion (i.e. when a = 2, [3 = 0), 
these boundary conditions are of a nonlocal nature; as 
a consequence, an analytical investigation of the mean 
first passage problem becomes very demanding and cum- 
bersome. In this work, we therefore restrict ourselves 
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predominantly (for a < 2) to detailed, precise numerical 
simulations. 

The setup for the studies is schematically depicted in 
Fig. n The stochastic motion of a free Levy particle is 
confined in an interval specified by the two boundaries 
Bi,i?2- The dynamics of the LFs trajectories derived 
from Eq. J2J) are consequently confined to this state space 
in-between B\ and Bi- 



x=0 



x=L 
B 2 



FIG. 1: Setup for the investigation of the first passage time 
analysis of Levy white noise driven free Brownian motion be- 
ing confined between the two boundaries located at B\ and 
B2, respectively. 



0.25 




0.05 




FIG. 2: Left panel: Typical sample trajectories of confined, 
normal diffusion, i.e. a — 2, with two absorbing AA bound- 
aries. Right panel: Confined normal Brownian motion be- 
tween a left-sided reflecting boundary (R) and a right-sided 
absorbing boundary (A), denoted as RA in the text. The 
particles start out at midpoint x/L — 0.5, undergoing state- 
continuous stochastic motion with evolving time t. 



A. The test case: normal diffusion 

For a = 2 (see Fig.|2J), white Levy stable noise becomes 
equivalent to Gaussian white noise and the corresponding 
Langevin eqution @ describes a free Brownian motion 



(Wiener process) for which the probability density func- 
tions (PDFs) of the first passage times distributions are 
known explicitly from the literature 0, Q . Here, due to 
reasons of convenience, instead of estimating the PDFs 
itself we simulate the equivalent cumulative distribution 
functions (CDFs) of first passage times. For the sake of 
simplicity it has been assumed that the left boundary 
is located at x = and the other one at x = L. In 
the following we typically use, although not exclusively, 
for the initial condition, xq, the center of the interval, 
xq = L/2. In the case when both B\ and B2 are absorb- 
ing boundaries, denoted as AA, the cumulative distribu- 
tion function J- of first passage times has been obtained 
by integration of the first passage time density f(t), i.e. 

= fo WW Hi US: 



T{t) = 1- 



2 " 1 

TT ^ 

J = l 



cos j7r . j7rao 
sm exp 

.7 L 



JTW_ 

L 



(6) 

Analogously, for a left absorbing boundary (A) and a 
reflecting right boundary (R) , denoted as AR, the corre- 
sponding CDF of first passage times reads 0, 0, IT(i| 



* w -1 

Ht) - -Y — 



1 — exp 



^(2J + I)*^(gj + W-*Q) 



2L 



(2j + l)an 
2L 



(7) 



Finally, for a left reflecting boundary (R) and a absorbing 
right boundary (A), denoted as RA, the corresponding 
CDF of first passage times reads 0, 0, ITo| 



7T ^ 2j 

j=0 J 



1 — exp 



. (2j + l)7T {2j + l)TTX 

■ sm cos ■ 



2L 



2L 



(8) 



where the scaling parameter a stands for the amplitude of 
the noise intensity for the additive white Gaussian noise; 
i.e. {((t)((s)) — 2a 2 5(t — s). Note, that for a particle 
starting its motion at xq = L/2, J-(t) for a RA boundary 
setup is the same as for the symmetrically chosen AR 
situation, cf. the formulae above. To test our employed 
software, we evaluated numerically the first passage times 
distributions with the same boundary conditions as spec- 
ified above. Eq. (0) has been integrated numerically for 
Bi = 0, B 2 = 5, x = 3.5 with At = ICT 5 and aver- 
aged over N = 3 x I0 4 realizations. Figure |3| depicts 
the numerical results along with their corresponding an- 
alytical expressions for the corresponding compound dis- 
tribution functions (CDFs) of first passage times for a 
Levy-Gaussian case (a = 2). We find a perfect agree- 
ment between the theoretical and the numerical results. 
On purpose, we have chosen here an asymmetric starting 
point in order to impose an explicit difference between 
the AR and the corresponding RA setup, respectively. 
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FIG. 3: First passage time statistics for free normal diffusion 
(a = 2) occurring between two boundaries located at x = 
and x — 5. The cumulative first passage time distribution 
T{t) is depicted versus the first passage time variable t. The 
initial motion starts out at xo = 3.5 and the noise strength is 
a 2 — 5. The chosen time step is At = 10 -5 and the results 
have been averaged over N — 3 x f 4 realizations. The simula- 
tion results are presented for AA, RA and AR configurations 
(from the top to the bottom) along with lines representing 
the analytical formulae for these cases. The numerical data 
match perfectly the analytical results, plotted on superim- 
posed lines. The inset depicts this agreement on an extended 
time span up to time t — 20. 



B. Confined Levy— Brownian motion 

After having tested the numerical algorithm, we next 
study the mean first passage time for confined Levy- 
Brownian motion on a finite interval. 



MFPT for symmetric Levy noise 

Using the discussed simulation procedure we start out 
with the case of confined Levy-Brownian motion with 
symmetric stable Levy noise, i.e. we set (3 — 0. Exem- 
plary trajectories are depicted in Fig.0J for the stability 
index a = 1.4. In clear contrast to normal Brownian mo- 
tion we detect discontinuous jumps, characterizing ran- 
dom jumps over large distances (note the horizontal ex- 
cursions in Fig. 0}. As a consequence, the boundaries 
at i?i,i?2 are typically not hit, but rather crossed in 
a flight-like manner, being characteristic for Levy dis- 
tributed jumps. In fact, for a < 2 large excursions of the 
trajectory are more probable than for normal Brownian 
motion. In effect, a test "particle" can skip over the bor- 
der, i.e. it can escape from the domain of motion via a 
single jump. 

As discussed in the literature 0, y, Hj, the scaling 
nature of the jump length PDFs causes a clustering of 
Levy flights. Random localized motion is occasionally 
interrupted by long sojourns on all length scales and, ad- 
ditionally, there are clusters of local motion within clus- 





FIG. 4: Sample trajectories for AA (left panel) and RA (right 
panel) with the two boundaries located at Bi = and B2 — L 
for symmetric Levy flights with a = 1.4,/? = 0. For a < 2 
large excursions of the trajectory are more probable than for 
the normal Brownian motion. A test "particle" can skip over 
the border, i.e. it can escape from the domain of motion with 
a single jump. 



ters. Anomalous trajectories of Levy flights with stability 
index a < 2 influence also the boundary condition for the 
mean first passage time and its corresponding statistics. 

From the perspective of a random walk approxima- 
tion to the Langevin equation in the behavior of a 
Levy walker is drastically different for Levy jump length 
statistics as compared with a traditional Gaussian case: 
The increments of the normal (state-continuous) diffu- 
sion process x{t) are characterized by statistics which 
excludes with very large probability the occurrence of 
long jumps. Therefore, the walker is more probable to 
approach - and eventually hit - a point-like boundary 
(cf. Fig. |2J). In contrast (see Fig.^J, with the Levy jump 
statistics (cv < 2), a meandering particle may easily cross 
the local boundary during its long jump and may recross 
into the finite interval many times, unless the particle 
is immediately absorbed upon crossing for the first time 
the boundary B\ or B2, respectively. This brings about 
a formulation of the boundary condition that necessarily 
must be non-local in nature. 

Indeed, the first passage time problem for Eq. (J2J can 
also be rephrased in terms of the fractional Fokkcr-Planck 
equation [H |M HI IH |H E3- For a free Levy- 
Brownian motion this equation with = assumes the 
form: 



dpjx, t) 
dt 



d\i 



■p(x,t) 



(9) 



The boundary conditions for the first passage time prob- 
lem associated with two absorbing boundaries at B\ and 
B 2 are now nonlocal; reading, 

p(x, t) = for x < B 1 and p(x, t) = for x ^ B 2 . (10) 

Due to discontinuous character of trajectories of Levy 



5 



processes we note that the usual form of boundary con- 
ditions, i.e. p(x — 0, t) — and p(x = L,t) = 0, in- 
correctly employed in the literature jH, 13, l3?iL l3ll | is 
expected to lead to erroneous results; such a boundary 
condition does not account for the fact that the process 
can skip the location of the boundary without hitting it. 
The corresponding MFPT can be numerically integrated, 
yielding 



t dt p(x, t)dx, 



MFPT 



which after a partial integration equals 



MFPT = I dt p(x,t)dx 

10 JBx 



(11) 



(12) 



This expression has been tested versus the statistical 
definition based on histogram analysis of the correspond- 
ing trajectories Eq. (J2J, for which the appropriate bound- 
ary conditions were imposed. 

Likewise, we also have studied the case where one of 
the two boundaries becomes reflecting. The reflecting 
boundary is imposed on the trajectory simulations via 
an infinite high hard wall, yielding immediate reflection 
|32| . For the numerical implementation of this reflecting 
case see below Eq. [5] In Fig. [5] we compare our numerical 
results for the case of two absorbing boundaries, case AA, 
and as well, for the symmetric situation, AR = RA, of a 
reflecting boundary and an absorbing boundary with the 
initial starting value chosen at midpoint. 
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FIG. 5: Mean first passage time versus the stability index a of 
confined Brownian motion driven by stable symmetric Levy 
noise. We depict the case of two absorbing boundaries, AA 
(left panel), and the situation of a reflecting and absorbing 
barrier, AR/RA (right panel). The initial starting point has 
been chosen at midpoint xo = 2, yielding identical results for 
AR and RA, respectively. The simulations parameters are: 
(3 — 0, a — \/2, time step At — 10~ 4 , number of realizations 
TV = 2 x 10 2 . The boundaries arc located at Bi = and 
B 2 = 4. The theory-result from Ref. ^ (Eq. (39) therein) 
is plotted as a dashed line. 



Fig. [5] depicts the numerical results for the MFPT for 
symmetric Levy stable noise after implementing numer- 
ically the appropriate boundary conditions as discussed 
above. In the left panel of Fig. the numerical results for 
the MFPT are compared with the theoretical findings for 
the MFPT = T in this superdiffusive case, reading from 
Eq. (39) in Refs. HEl, 



T 



ttD 



a oo 



(-1)' 



^ (2m + l) 1+ «' 

rn— v ' 



(13) 



where D a = a a . 

The maximal MFPT is assumed for a ~ 1. The behav- 
ior of the MFPT assumes a bell-shaped behavior around 
a = 1 ; it reflects the different interplay of the probability 
of finding long jumps versus a decreasing stability index 
a < 2, implying a decreasing noise intensity a a . Put dif- 
ferently, the occurrence of long jumps beyond the finite 
interval is dominating the escape over a decreasing noise 
intensity for for a < 1. 

The expected discrepancy between these analytical re- 
sults and numerical estimation is due to assumed local 
boundary condition of vanishing probabilities p{x = 0) = 
p(x = L) = 0, which are correct only for normal Brown- 
ian motion; i.e. a — 2. Indeed, our numerical result just 
coincides precisely at this very special value. In our simu- 
lation of the Langevin equation in (j2J the whole exterior 
of the prescribed interval (E>i = 0, P>2 — 4) is absorb- 
ing throughout, while for the analytical calculations the 
flawed, point-like boundary conditions are assumed |22|. 

Clearly, the extension of the absorbing regime to the 
whole two semi-lines outside the confining interval yields 
the physically correct value that accounts for the escape 
of the particle from the interval via long jumps. The 
numerical results for a < 2 systematically exceed the 
theory-result of Ref. [23 • 

The right panel of Fig. [S] depicts the case for the MF- 
PTs with AR and RA boundaries, respectively. Due 
to symmetric chosen initial condition and the symmet- 
ric stochastic driving the results for AR and RA become 
identical. 



MFPT for asymmetric Levy noise 

Asymmetric Levy noise is characterized by a non- 
vanishing skewness parameter f3 ^ (see the exemplary 
probability density functions in Fig. EJ. Fig. presents 
the results of MFPT evaluations for fully asymmetric 
Levy stable noise driving with \(3\ = 1. The left panel dis- 
plays results for the Levy-Brownian motion with (3 = 1 
within the AA boundaries setups. Due to the imposed 
symmetry of a starting point, the results coincide with 
those for the Levy noise driving with (3 = — 1. The right 
panel depicts the results for AR boundaries with (3 = ±1. 
As can be inferred by inspection of Fig. due to the 
skewness character of the Levy stable distribution, the 
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results for AA boundaries with a < 1 and (3 = 1 are the 
same as for AR boundaries with a < 1 and f3 = — 1. The 
effect is caused by a visible shift of the probability mass 
apart from x — to the left (or to the right) for skewed 
distributions with a < 1. In contrast, these results differ 
for a > 1 where both boundary setups, i.e. AA and AR, 
lead to different values for MFPTs. 
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FIG. 6: Probability density functions for the a-stable vari- 
ables x with a — 1.5 (left panel) and a — 0.9 (right panel). 
The symmetric case is for (3 = 0, while [3 — ±1 corresponds 
to asymmetric densities. Note the differences in the positions 
of the maxima for a < 1 and a > 1. Most importantly, the 
support of the densities for the fully asymmetric cases with 
(3 = ±1 and a < 1 (right panel) is covering not the whole 
axis; it assumes only negative values for /3 = — 1 and only 
positive values for (3—1. This in turn causes the discontinu- 
ous behavior depicted with Fig. |7| 



passage time distributions (CDFs) of the first passage 
times as governed by Eq. @ for the AA boundaries with 
a symmetric starting point, i.e. xq = L/2, are invariant 
under the transformation (3 — * — /3. Therefore, the quan- 
tities derived from the numerically determined CDFs are 
symmetric around f3 = 0, cf. the left panel of Fig- EI Fur- 
thermore, because of the chosen symmetric initial condi- 
tions the cumulative first passage time distributions eval- 
uated for the RA boundaries setup can be constructed 
from the appropriate distribution for the AR case by ex- 
changing j3 with — (3. It is caused by the fact that the 
asymmetry induced by a nonzero skewness parameter can 
be compensated by an exchange of boundaries, cf. left 
panel of Fig. [SJand Fig. 
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FIG. 8: MFPTs versus the skewness parameter (3 for AA (left 
panel) and AR (right panel) boundaries for Levy-Brownian 
motion with a — 0.5. The simulations parameters are like in 
Fig.El 
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FIG. 7: MFPTs versus the stability index a for asymmetric 
Levy noise for different boundary conditions: A A (left panel), 
AR (right panel). The simulations parameters are as in Fig. El 
Note that the results for the AR setting with (3 = —1 and 
a < 1 are the same as for the AA settings with j3 ± 1 and 
a < 1. 

As can be intuitively expected, the cumulative first 



C. First passage time statistics 

The simulation data for the mean first passage time 
yield as well the results for the cumulative first passage 
time distribution functions (CDFs). Typical such dis- 
tributions are depicted in Fig. [5] for symmetric a-stable 
noise and in Fig. 1101 for asymmetric stable Levy white 
noises with a stability index a — 1.5. With symmet- 
ric noise we observe that the escape, starting out at 
midpoint, is speeded up with two absorbing boundaries 
AA as compared to the situation with AR=RA. The re- 
flecting boundary clearly slows down the ultimate escape 
from the finite interval. 

For escape driven by asymmetric white Levy noise, 
see Fig. El we note that for the case of two absorb- 
ing boundaries AA (left panel) the CDFs become iden- 
tical for (3 = —f3, given the midpoint starting value. We 
also can detect a more rapid saturation with (3 ^ as 
compared to the fully symmetric situation with f3 = 0; 
implying a somewhat faster escape scenario. This fact 
originates from the skewness in the distribution of jump 
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values, implying a faster escape towards the correspond- 
ing absorbing boundary. 

With an AR boundary setup the situation becomes 
more intricate. Now, depending on the choice of the 
asymmetry parameter the escape can be enhanced, re- 
flecting the skewness of the corresponding stable distri- 
bution and the relative character of the boundary setup. 




t t 



FIG. 9: Cumulative first passage time distribution function 
for symmetric Levy-Brownian motion on a restricted interval, 
i.e. f3 = 0, with a = 1.9, 1.0, 0.5 for AA (left panel) and AR 
(right panel) boundaries. The results have been simulated for 
a = \/2, At = 10~ 4 and have been averaged over N = 2 x 10 2 
realizations. The boundaries Bi , B2 are located at x = and 
x = 4. The particle starts out at midpoint; i.e. xq = 2. 



values only (3 = — 1. Therefore, the results for a < 1 
and \(3\ = 1 for AA boundaries are the same as for AR 
((3 = —1) or RA ((3 — 1) boundaries. This effect is clearly 
visible for the Levy-Smirnoff case depicted in Fig. ^2 
Now, the results for AA and RA boundaries are the same, 
as expected. A similar behavior is observed in Fig. 
above, where the values of the MFPTs for a < 1 for AA 
and AR boundaries agree. 
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FIG. 11: The cumulative first passage time distribution J-(t) 
versus first passage time t for Levy-Smirnoff driven Levy- 
Brownian motion, i.e. a = 0.5 and /3 = 1 and a scale param- 
eter a — \/2- Because the Levy-Smirnoff statistics assumes 
only positive values the distributions J-(t) for AA and RA 
configurations become the same. The simulations parameters 
are as in Fig. |S] 




t t 



FIG. 10: Cumulative first passage time distribution function 
versus first passage time t for Levy noises with a = 1.5 and 
P = -1.0,0.0,1.0 for AA (left panel) and AR (right panel) 
boundaries. The simulations details are the same as in Fig. El 
Note that the results for the AA settings with /3 = ±1 overlap. 

This complexity is elucidated further in panel 1111 for 
the case of stable Levy-Smirnoff noise where a = 0.5, (3 — 
1. The CDF T{t) is depicted in Fig. [TT] For a < 1 
and |/3 1 = 1 Levy stable distributions are totally skewed, 
taking on either only positive values (/3 = 1) or negative 



Survival probability: comparison with Sparre- Andersen 
scaling 

The survival probability <S(t) = 1 — T{t) for symmetric 
Levy-Brownian motion on the finite interval is displayed 
in Fig. 1121 The motion on a confined support leads to 
an exponential decay of the survival probability with the 
steepness of the slope depending on the stability index a. 
Notable differences can be observed for various bound- 
ary setups: For motion occurring between two absorbing 
boundaries, the survival probability decays faster than 
for the case when one of the boundaries is taken as re- 
flecting. It is also useful to emphasize that deviations 
from the exponential behavior of the survival probability 
can be observed in systems subjected to both, dichotomic 
and Levy stable noises |33l l34| , respectively. 

Finally, for the systems driven by Levy stable noises 
we also tested the Sparre- Andersen scaling behavior on 
an infinite half line. According to the Sparre-Andersen 
theorem for a free stochastic processes driven by sym- 
metric white noises, the first passage time densities 
process from the real half line asymptotically behave like 
i -3 / 2 . Consequently the survival probability, i.e. the 
probability of finding a particle starting its motion at 
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FIG. 12: Survival probability S(t) = 1 - F{t) for confined, 
symmetric Levy-Brownian motion (5 = on the interval 
Bi = 0, B2 — 4 with absorbing half- line for x < and ab- 
sorbing (or reflecting) half-line for x > for various values of 
the stability index a and midpoint initial conditions at xo = 2. 
The left panel depicts the results for AA boundary setups; the 
AR cases are displayed with the right panel. The solid line in 
both data sets represents the power law t~ ll/2 which describes 
an asymptotic behavior foreseen by the Sparre- Andersen the- 
orem 1131 for diffusion on semi-infinite intervals. 



xo > in the real half line, scales like i -1 / 2 . In Fig. 1131 
the survival probability S(t) = 1 — Tit) is depicted for 
various stability indices a and various initial conditions 
xq. It is clearly visible that the survival probability S(t) 
behaves like a power law with the exponent (—1/2), as 
predicted by the Sparre- Andersen theorem. For the test- 
ing of the Sparre-Andersen theorem the whole negative 
half line was assumed to be absorbing. 
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FIG. 13: Survival probability S(t) = 1 - F(t) for free Levy 
flights on the half line with absorbing boundary for x < 
for various values of stability index a and various initial con- 
ditions Xq. Remaining simulation parameters like in Fig. |5] 
The survival probability nicely fits i~ 1//2 slope predicted by 
Sparre-Andersen scaling. 



IV. CONCLUSIONS 

With this work we investigated the problem of the 
mean first passage time and the first passage time statis- 
tics for Markovian, Levy-Brownian motion proceeding 
on a finite interval. The intricate problem of setting up 
the proper boundary conditions for absorption and re- 
flection are discussed with possible pitfalls being pointed 
out. In particular, it has been demonstrated by numer- 
ical studies that the use of the commonly known, local 
boundary condition of vanishing flux (in case of reflec- 
tion) and vanishing probability (in case of absorbtion), 
valid for normal Brownian motion (i.e. a = 2) no longer 
apply for Levy white noise. This is so, because the large, 
jump like excursions of Levy flight increments causes non- 
continuous, i.e. discontinuous sample trajectories. This 
in turn requires the use of nonlocal boundary conditions. 
It is presently not known how these nonlocal boundary 
conditions can be recast as in equivalent form as a mod- 
ified, locally defined differential condition involving the 
statistical quantities of interest (the MFPT and the first 
passa ge t ime densities) at the location of the boundary 
alone JM H^. 

For symmetric, a stable Levy white noise driven Brow- 
nian motion we find a bell-shaped, non-monotonic behav- 
ior of the MFPTs for absorbing-absorbing boundaries, 
with the maximum being assumed for a = 1. In contrast, 
for reflecting-absorbing boundaries we find a monotonic 
increase. With asymmetric Levy white noise, i.e. with 
a non-vanishing skewness parameter (3^0 the MFPT 
results in an even more complex behavior. As a function 
of the stability index a the MFPT can exhibit a discon- 
tinuous behavior, see the right panel in Fig. [7] 

In addition we have studies also the statistics of the 
first passage times in terms of the cumulative first pas- 
sage time distribution function (CDF) and the corre- 
sponding survival probability. In this context, we also 
belabor the role of a finite support and different bound- 
ary setups for Levy white noise driven Brownian motion 
for the universal scaling law of Sparre-Andersen. While 
the restricted Levy-Brownian motion exhibits an expo- 
nential behavior on finite intervals, the crossover to the 
universal i -1 / 2 - power law, being valid on the half-line for 
all stability indexes is assumed for large intervals only. 



APPENDIX A: a-STABLE RANDOM VARIABLES 

The random variables £ corresponding to the charac- 
teristic functions and ffl can be generated using the 
Janicki-Weron algorithm [la . Il7| . For a^l, their rep- 
resentation reads 



sin(a(V + C a , fi )) 



(cos(7)) = 

cos(V - a(V + C a ,p)) 
W 



(Al) 
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with the constants C, D given by 

arctan(/3tan(^)) 



C a ,p — 



a 



(A2) 



In contrast, for a 
distribution 



l,j3 = one obtains the Cauchy 



- £ 'i,o(C; cr ,M) = 



7r(C-/^) 2 + a 2 



(A6) 



cos ^arctan 



( n , 7TQ; N 

(0tan( — ) 



(A3) 



For a = 1, the random variable £ can be calculated from 
the formula 



2rr 
n 



(- + f3V)t & n(V) 



%Wcos(V) 



(A4) 



In the above equations V and W denote independent 
random variables; namely, V is uniformly distributed in 
the interval (— f , f) while VF is exponentially distributed 
with a unit mean [16| . The numerical integration scheme 
has been performed for /i = with the increments of Ax 
(see Eqs. (0) and 10) sampled from the strictly stable 
distributions 

The analytical expressions for stable probability dis- 
tributions L a .p((;<7,fi) are known in few cases only: For 
a = 0.5, f3 = 1 the resulting distribution is the Levy- 
Smirnoff one; i.e., 



£i/2,i(C;ct,/i) = (^) 2 ((^^y 



x exp 



2(C - M) 



(A5) 



The familiar case with a = 2 with arbitrary (3 yields the 
Gaussian PDF. The prominent characteristic feature of 
the distributions L a ^(<^; a, fi) is its existence of moments 
up to the order a, i.e. the integral /_ -L ai( a(£; a, /x)£ a d£ 
is finite. This statement results in the conclusion that the 
only stable distribution possessing a finite second moment 
is the Gaussian; for all other values of a the variance of a 
stable distribution diverges, and for a < 1 even the first 
moment does not exist. 
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